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1 Introduction 

Recently the considerable advance in the rigorous derivation of quantum kinetic equations in 
the mean field scaling limit [1], such as the nonlinear Schrodinger equation and the Gross- 
Pitaevskii equation pHTT] is observed. In the scaling limit approach to the construction of 
kinetic equations one of the essential assumptions is that initial data has to satisfy a chaos 
property, i.e. the initial state is completely determined by a sequence of products of the one- 
particle density operators [1], [12]. At the same time it is well known that, for instance, the 
equilibrium state of the Bose condensate is characterized by correlations of particles in contrast 
to the gaseous state [13]. Thus, the validation of quantum kinetic equations in case of presence 
of correlations at initial time, in particular for condensed states, remains an actual problem so 
far. 

In the paper we consider the problem of potentialities inherent in the description of the 
evolution of states of many-particle systems in terms of a one-particle density operator in case 
of presence of correlations of particles at initial time. We prove that, if initial data is completely 
specified by a one-particle marginal density operator and initial correlations, then all possible 
states of infinite-particle systems at arbitrary moment of time can be described without any 
approximations within the framework of a one-particle density operator governed by the kinetic 
equation, which coefiicients are determined by the initial correlations. Then the actual kinetic 
equations that describe the evolution of interacting particles in the condensed states, are derived 
in appropriate scaling limits on the basis of the established quantum kinetic equation. 

We consider a quantum system of a non-fixed (nonequilibrium grand canonical ensemble) 
number of identical (spinless) particles obeying Maxwell-Boltzmann statistics in the space M^. 
We will use units where h = 2'Kh = 1 is a Planck constant, and m = 1 is the mass of particles. 

Let l-ihe a, one-particle Hilbert space, then the n-particle space is a tensor product of 
n Hilbert spaces 1-L. The Hamiltonian of n-particle system is a self-adjoint operator with 
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domain V{Hn) C Tin- 

n n 

H^ = Y,K{i)+e J2 H^u^2), (1) 

i=l n<j2=l 

where K{i) is the operator of a kinetic energy of the ith particle, ^2) is the operator of a 
two-body interaction potential and e > is a scaling parameter. The operator K{i) acts on 
functions ■0n, that belong to the subspace Lq(]R^") C T>{Hn) C L^(]R^") of infinitely differen- 
tiable functions with compact supports according to the formula: K{i)'i{jn = —^Ag-ipn- Cor- 
respondingly we have: $(zi,Z2)'0n = ^{ih, Qi2)'^n, and we assume that the function ^{qii^Qi^) 
is translation-invariant bounded function which is symmetric with respect to permutations of 
arguments. 

Let £^(7{„) be the space of trace class operators /„ = /n(l, . . . , n) G that satisfy 

the symmetry condition: . . . , n) = /n(^i, ■ ■ ■ ,in) for arbitrary (ii, . . . , z„) G (1, . . . , n), and 

equipped with the norm: ||/ni|£i(w„) = Tri_ . . . , n)|, where Tri_ „ are partial traces 

over 1, . . . ,n particles. We denote by £q('H„) the everywhere dense set of finite sequences of 
degenerate operators with infinitely differentiable kernels with compact supports. Let -C('H„) 
be the space of bounded operators defined on the Hilbert space T-Ln- 

On the space Sl}{J^y) the one-parameter family of operators 

6^n(-t)/n = e-^*^"/„e**^". (2) 

is defined and it is an isometric strongly continuous group which preserves positivity and self- 
adjointness of operators. For / G ^KTin) there exists the infinitesimal generator of this group 

lim - /„) = -i{Hnfn - fnHn) = -J^nfn, (3) 

where if„ is Hamiltonian ([T]) and the operator: —i{Hnfn — fnHn) is defined on the domain 

V{Hn) C Kn. 

Let us denote by {Y} the set consisting of one element Y = (l,...,s). We define the 
(1 + n)f/i-order {n > 0) cumulant of groups of operators ([2]) as follows 

^,U-t,{Y},X\Y)= Yl (-1)"''"'(|P|-1)! n ^l^(^OlH'W), (4) 

P:({Y},X\Y)=[J^Xi X,CP 

where {Y} is the set consisting of one element F = (1, . . . , s), is the sum over all possible 
partitions P of the set {{Y}, X \Y) = {{Y}, s + 1, . . . , s + n) into |P| nonempty mutually 
disjoint subsets Xi C {{Y}, X \ Y) and 9 : {{Y}, X \ Y) ^ X is the declasterization mapping. 

We will consider initial state which is given by the following sequence of marginal density 
operators 

n 

F{t)U=o = (1, F^{1),92{1, 2)F?{1)F?{2), ^„(1, n) J] F^{^), • • • ) ' 

i=l 

where the bounded operators gn G £(7^n), > 2, are specified initial correlations. Such initial 
data is typical for the condensed states of quantum gases, for example, the equihbrium state of 
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the Bose condensate satisfies the weakening of correlation condition with the correlations which 
characterize the condensed state [I3]. We remark that nonequilibrium dynamics of correlations 
was constructed in [H]. 

The evolution of all possible states is described by the sequence of marginal density operators 
F{t) = (l, 1), . . . , Fs{t, 1, . . . , s), . . . ) governed by the initial- value problem of the quantum 
BBGKY hierarchy 

i=i ii<j2=i 

s 

+e J2 Tr,+i{-MntU, s + Y,s + 1), (5) 

s 

F,{t,Y)\t=o=9s{Y)l[F^{t), s>l, (6) 

i=l 

where on /„ G ^liT-Ln) C S^^iTin) the operators {—J\f{j)) and {—N'mtUi^h)) ^-^e defined by the 
corresponding formula 

(-A/'(j))/n = -^{K{j)fn - fnK{j)), (7) 
(-A/'int)(jl, j2)/„ = -i($0l,j2)/„ - /n$(jl,i2)). 

If f;i(^) < and max„>2||5'„||£(-^„) < oo, then for t G M a nonperturbative solution of 
the Cauchy problem ([5])- ([6]) is given by the expansion [15j 

oo ^ s+n 

F,{t, Y) = Y,^_ Tr,+i,...,,+„ 2ti+„(-t, {Y},X\ Y)g,+,,{{Y}, X \ y) J] F^{t), (8) 

n=0 ■ i=l 

where 2li+„(— t) is the (1 + n)t/i-order (n > 0) cumulant (jlj of groups of operators ([2]) and 
gi+n{{Y}, X\Y) is the correlation operator of clusters of particles [T4], in particular case n = 
it has the form: gi+o{t, {Y}) = Ep:y=u,x,nx,cP^l^.l(^' ^')- 

In the paper it was proved that, if initial data is completely determined by a one-particle 
marginal density operator, i.e. gi+n = /, n > 1, in initial data ([6]), then states given in terms 
of the sequence of marginal density operators (IH]) can be described within the framework of the 
sequence F{t \ Fi{t)) = (1, Fi{t),F2{t \ Fi{t)), . . . , | . . .) of explicitly defined func- 

tionals Fs{t \ Fi{t)), s > 2, of the solution Fi{t) of the generalized quantum kinetic equation. 

We extend this statement to case of initial data ([6]) and construct the mean field (self- 
consistent field) asymptotics of a solution of the Cauchy problem of the corresponding gener- 
alized quantum kinetic equation. 



2 The generalized quantum kinetic equation in case of 
correlated initial state 

We reformulate the Cauchy problem ([5])- ([6]) as the new Cauchy problem for a one-particle 
density operator governed by the generalized quantum kinetic equation and the sequence of 
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explicitly defined marginal functionals of the state Fs[t, Y \ i^i(t)) , s > 2, which are determined 
by the solution Fi{t) of such Cauchy problem. With this aim we introduce the following kinetic 
cluster expansions of cumulants @ of groups of operators ([2]) 



s+n 



2ti+„(-t, {Y}, X \ Y)g,+4{Y}, X \ y) J] t) = (9) 
n! 

1 

{n — ni)\ 



i=l 

n . 

ETl' 
7 ^6i+n-ni(t, {Y}, s + l,...,s + n-ni) X 

ni=0 

s+n— ni |Ds+n| ^ 

Ds+n 

X n 9i+\Xk\i^k,Xk)^i{t,ik)Y[^iitJ), 

k=i jez 

where X\Y = {s + 1, . . . ,s + n) and J2Ds+n-z=[j Xi '^^^^ possible dissections J^s+n 

of the linearly ordered set Z = {s + n — ni + 1, . . . , s + n) on no more than s + n — ni linearly 
ordered subsets. We give a few examples of recurrence relations ([9]) 

^,{t,{Y}) = (3^{t,{Y}), 

s 

{Y},s + 1) = &2it, {Y}, s + l) + e^{t, {Y}) Mt, H,s + l), 

h=i 

2t3(t, {Y},s + l,s + 2) = 03(t, {Y}, s + 1, s + 2) + 

s+l s 

+2!(52(t, {Y},s + 1)Y ^2{t, ti,s + 2) + (3,{t, {Y}) ( ^ M^, ii,s + l,s + 2) + 



ii = l ii = l 



+2! Y Mt,H,s + l)^2it,i2,s + 2)), 



l=ii<i2 

where the operator 2li+„(t) is the (1 + n)t/i-order scattering cumulant 



s+n 



{Y}, X\Y) = {Y}, X\Y) J] 2ti(f, t), (10) 

1=1 

In terms of scattering cumulants ( ITOl) the solutions n > 0, of recurrence 

relations ([9]) are represented by the following expansions 

n n n-ni-...—nk-i 

«.UMn^\n = n!5:(-i)'5:... ^^-^x (n) 

fc=0 n.i = l nfc=l ^ 

X^l+n-ni-...-n^it, {Y} , S + 1, . . . , S + U - Ui - . . . - Hk) X 
k ^ s+n-ni-...-nj ^ 

D, I < s + n — ni — • ■ ■ — n,- 
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where XId z =11 x ^^^"^ possible dissections of the hnearly ordered set Zj = 

{s + n — rii — . . . — rij + 1, . . . , s + n — TLi — . . . — fij^i) on no more than s + n — ni — . . . — rij 
hnearly ordered subsets. For example, 

s 

(5i(t, {Y}) = 2ti(t, {Y}) = {Y})gr{{Y}) J] 2li(t, z), 

i=l 

s 

(32{t, {Y}, s + l)= Ql2{t, {Y},s + 1)- 2li(t, {Y}) Mt, ^,s + l), 

i=l 

iS^it, {Y}, s + l,s + 2)= 2t3(t, {Y}, s + l,s + 2)- 

s+l s 

-2!2l2(t, {Y}, s + l)J2Mt,^uS + 2)- 2ti(t, {Y}) ( ^ ^,{t, t,,s + l,s + 2)- 

il=l ii=l 
s s+l s 

-2!^^2i2(t,ii,s + l)2l2(t,Z2,s + 2) + 2! J] 2t2(t, ii, s + l)2l2(t, ^2, s + 2)). 

«i=i*2=i i=n<i2 

As a result of the application of expansion ([9]) to solution ([8]) we establish that the one-particle 
density operator Fi{t) is governed by the following generalized quantum kinetic equation 

j^F,{t,i) = -mmit,i)+ (12) 

00 ^ n+2 

+e Tr2(-Ari,t(l, 2)) ^ ^Tr3„„,„+2 ^i+n(t, {1, 2}, 3, . . . , n + 2) J] F,{t, i), 

n=0 ■ i=l 

Fmt=o = F^- (13) 

where the operator (— A/'int(li 2)) is defined by formula ([7]), and the (1 + n)t/i-order generated 
evolution operators 0i+„(t), n> 0, are determined by expansion ( ITTj) over scattering cumulants 
( !T0|) . In equation f|T2|) the collision integral series converges under the condition that [16] : 
||Fi(t)||.i(^) <e-^ 

The global in time solution of initial- value problem f ll2p -f ll3p is determined by the following 
expansion 



n+l 



1) = 5Z ^ Tr2,...,i+„ 2li+„(-t, 1, 2, . . . , n + l)(7i+„(l, 2, . . . , n + 1) J] ^^(z), (14) 

n=0 ■ 1=1 

where 2li4.„(— t) is the (1 + njth-ordei cumulant (jl]) of groups of operators ([2]). The series ( IT^ 
converges under the condition that [16]: ||F°||£i(^) < e~^°(l + e^^)^^. 

Correspondingly having applied expansion ([9]) to solutions ([H]) in case of s > 2, in terms of 
solution expansion ( !T^ . they are represented as the following marginal functionals of the state 

00 s+n 

Fs{t, Y I F,{t)) = J2^_ Tr,+i,„„,+„ {Y},X\Y)l[ F,{t, z), (15) 

n=0 ■ i=l 

where the (1 + n)t/i-order generated evolution operator 0i+„(f) is determined by expansion 
(fTTj) over scattering cumulants (fTOj) . The series (fT5l) converges under the condition that: 
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\\Fi{t)\\zi(H) < e-(^'+2). The constructed marginal functionals of the state ( ITSll characterize 
the correlations of states of quantum many-particle systems. 

We note that the average values of the nonadditive-type marginal observables are determined 
in terms of marginal functionals of the state ( !T5|) and in case of the additive-type marginal 
observables they are determined by solution (fT^ of the generalized quantum kinetic equation 

m- 

Thus, in case of initial data solution (|H]) of the Cauchy problem ([S])-® of the quantum 
BBGKY hierarchy and a solution of the Cauchy problem of the generalized kinetic equation 
f|T2|) - f|T3|) together with marginal functionals of the state f|T5|) give two equivalent approaches 
to the description of the evolution of quantum many-particle systems. The coefficients of 
generalized quantum kinetic equation f|T2|) and generated evolution operators (|TT1) of marginal 
functionals f llSp are determined by the operators of initial correlations. 



3 The mean field limit theorems 

In case of initial state involving correlations for an asymptotic perturbation of generated evo- 
lution operator (ITTl) in the mean field limit the following equality is valid 

lim||0i+„(t,{y},X\r)/,+„|| =0, n>l, 

and in case of the first-order generated evolution operator ffTTl) we have respectively 

s s 

hmJ\{0,{t,{Y})-]lg,{-t,z,)g,{{Y})]lg,{t,z,))f,\\^ 

il = l 42 = 1 

In view that under the condition that: t < to = (2 ||$||£c^2) ||e F°||£i(-^)) ^, the series for 
e Fi (t) converges, then for t < to the remainders of solution series (fT4|) can be made arbitrary 
small for sufficient large n = hq independently of e. Then, using stated above asymptotic 
perturbation formulas for each integer n every term of this series converges term by term to 
the limit operator fi{t) which is represented in the form of the following expansion 

f,{t,l) = J2 dil... / rft„Tr2,...,i+ne?i(-t + ti,l)(-Mnt(l,2)) X (16) 

n=0 { i 
2 n 
X Yl Gli-h + t2, jl) . . . Yl Gli-tn-1 + tn, Jn-l) X 
ii=i in-l=l 
n l+n 1+n 

X ^ (-A/;nt(^n,l + ^)) n Gli-tn,Jn)9l+n{l,...,n+l)l[f^{l). 

in = l jn = l «=1 



For bounded interaction potentials series (ITBI) is norm convergent on the space under 

the condition that: t < = (2 \\^\\s.(H2)\\fi\\2^iH)) 

Thus, if there exists the limit /° G £,^{71) of initial data (fT3|l . namely 

Slh^i-/?IUH) = 0' 
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then for finite time interval t G (— ^0,^0), where to = (2 ||'^'||£(-H2)II/i II^M^)) ^' there exists the 
mean field limit of solution expansion f ll4p of the generalized quantum kinetic equation f ll2p : 

lim||eFi(t)-A(t)||^,(^^ = 0, (17) 

where the operator /i {t) is represented by series fll6p and it is a solution of the Cauchy problem 
of the modified quantum Vlasov kinetic equation 

|A(t,l) = -A/-(l)A(t,l)+ (18) 

2 2 

+Tr2(-A/;nt(l, 2)) n Sii-t^ ^i)9ii{h 2}) n 2), 

ii = l 12=1 

/i(t)k=o = /?. (19) 

Since a solution of initial- value problem (IT^ - fll3l) of the generalized kinetic equation con- 
verges to a solution of initial- value problem f|T8|) - f|T9|) of the modified quantum Vlasov kinetic 
equation as (|T71) . for marginal functionals of the state f|T5|) we correspondingly establish 



lim||e^F,(t,y I F,{t)) - l[g,{-t,t,)g,i{Y}) H ^2) J] /i(^'^')Lhw.) = ^^0) 

ji=i 42=1 i=i 

This equality means the propagation of initial correlations in time in the mean field limit. 



4 On mean field quantum kinetic equations 

Let us consider the pure states, i.e. the operator fi(t) = {i/jt) {i^t \ is a one- dimensional projector 
onto a unit vector \ipt) G "H and its kernel has the following form: fi{t,q,q') = ip{t,q)'i/j(t,q'). 
Then, in case of initial data, satisfying a chaos property, it holds 

lim||eF°-|V^o)(^olL.(^)=0, 
and in this case statement fl20|) reads 

s s 

lim II e'F,{t I Fi(t)) - J] ^i(-t, ^l)^7l(m) H ^2)|^*)(^r^ h^^n.) = 0' 

n=i 12=1 

where the function {ipt) G is a solution of the Cauchy problem of the nonlinear Hartree 
equation for initial data \ipo). In particular, we remark that in case of a system of particles, 
interacting by the potential which kernel is the Dirac measure, it reduces to the cubic nonlinear 
Schrodinger equation 

^^^(^,9) = ~A,i;{t,q) + \ij{t,q)\^ij{t,q). 
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In case of correlated initial state ([6]) the sufficient equation for the description of the pure 
state evolution governed by modified quantum Vlasov kinetic equation ( fT8|) is the Gross- 
Pitaevskii-type kinetic equation 

z-i{j{t, q) = --A.ijit, q) + J dq'dq"b{t, q, q; q' , q")^{t, q")r{t, q)i:{t, q'), (21) 

where the coupling ratio function b{t, q, q; q', q") is the kernel of the scattering length operator: 

Observing that in the macroscopic scale of the variation of variables, groups of operators 
(I2]) of finitely many particles depend on microscopic time variable e~^t, where e > is a scale 
parameter, the dimensionless marginal functionals of the state are represented in the form: 
Fs(^e~^t I Fi{t)Y As a result of the limit processing e — )■ we establish the Markovian kinetic 
evolution fl2T|) with the corresponding coefficient b{e~^t). 

5 Conclusion 

We generalized the concept of quantum kinetic equations for the kinetic evolution involving 
correlations of particle states at initial time [12], for instance, correlations characterizing the 
condensed states [I3]. The mean field scaling asymptotics of a solution of the generalized 
quantum kinetic equation of many-particle systems in condensed states was analyzed. These 
results can be extended to quantum systems of bosons or fermions . 

It should be emphasized that the kinetic evolution is an inherent property of infinite-particle 
systems. In spite of the fact that in terms of a one-particle marginal density operator from 
the space of trace-class operators can be described a system with the finite average number 
of particles, the generalized quantum kinetic equation has been derived on the basis of the 
formalism of nonequilibrium grand canonical ensemble since its framework is adopted to the 
description of infinite-particle systems in suitable Banach spaces as well. 

We note that one more approach to the construction of the kinetic equations in the mean 
field limit in case of the presence of correlations at initial time is based on the description of 
the kinetic evolution in terms of the evolution of marginal observables [T71IIB]- 

We remark that developed approach is also related to the problem of a rigorous derivation 
of the non-Markovian kinetic-type equations from underlaying many-particle dynamics which 
make possible to describe the memory effects of particle. 
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